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Patterns of convection in internally heated, self gravitating rotating spher- 
ical fluid shells are investigated through numerical simulations. While turbu- 
lent states are of primary interest in planetary and stellar applications the 
present paper emphasizes more regular dynamical features at Rayleigh num- 
bers not far above threshold which are similar to those which might be ob- 
served in laboratory or space experiments. Amplitude vacillations and spatial 
modulations of convection columns are common features at moderate and 
large Prandtl numbers. In the low Prandtl number regime equatorially at- 
tached convection evolves differently with increasing Rayleigh number and 
exhibits an early transition into a chaotic state. Relationships of the dynam- 
ical features to coherent structures in fully turbulent convection states are 
emphasized. 

1 Introduction 

Thermal convection in a rotating spherical fluid shell represents one of the funda- 
mental systems on which a considerable part of the understanding of many observed 
geophysical, planetary and astrophysical phenomena is based. It is generally believed 
that the magnetic field exhibited by celestial bodies is generated in their interiors 
by convection driven by thermal and compositional buoyancy. Cloud patterns and 
the differential rotation seen at the surface of the major planets are likely to be 
connected with convection in the deep atmospheres of those planets (Busse, 1976, 
1994). Considerable attention has therefore been devoted to the study of dynamical 
properties of rotating spherical shells. For recent reviews we refer to the papers by 
Zhang and Busse (1998) and Busse (2002). 

Because a radially directed buoyancy force can not be easily realized in Earth- 
bound laboratories, experimenters had to resort to the use of centrifugal accelera- 
tions (Busse and Carrigan, 1976) or use spacecrafts for their experiments (Hart et al. 
1986). Fortunately, it has turned out that the component of gravity perpendicular 
to the axis of rotation is much more important than the component parallel to it. 



1 



The centrifugal force is thus a convenient substitute since through the reversal of the 
applied temperature gradient the same buoyancy can be obtained as for the inward 
directed component of gravity perpendicular to the axis. Moreover, the centrifugal 
force increases with the distance from the axis just as the case of the gravity field 
of a spherical body of constant density. The use of electric fields (Hart et al. 1986) 
or magnetic fields with ferrofluids (Rosenzweig et al. 1999) leads to rather different 
radial dependences of the effective gravity. It must be admitted, however, that the 
absence of an axial component of gravity in experiments employing the centrifugal 
force leads to a thermal wind type flow in the basic state. But this can be minimized 
when the Coriolis parameter is sufficiently large or when a thin spherical shell is used 
as for example in the experimental realization of convection shown in figure [U 

In this paper results from a numerical study will be reported that have been 
obtained with essentially the same computer code as employed in earlier studies 
focusing on turbulent convection (Tilgner and Busse, 1997) and on the generation 
of magnetic fields (Grote and Busse, 2001). The mathematical formulation and 
the description of the numerical method given in section 2 can thus be kept brief. 
In section 3 a short introduction to the results of linear theory for the onset of 
convection will be given. In section 4 properties of columnar convection at Prandtl 
numbers of the order unity will be described and some animations of the time 
dependence will be presented. Equatorially attached convection flows are the subject 
of section 5. Their properties at finite amplitudes are studied here for the first time. 
Because of the required high numerical resolution in time as well as in space the 
low Prandtl number regime of fully nonlinear convection has posed a significant 
computational challenge until now. As will be shown the continuity of the convection 
heat flux requires that equatorially attached convection must be joined by columnar 
convection as the Rayleigh number increases significantly above its critical value. A 
concluding discussion will be given in section 6. 

2 Mathematical Formulation and Numerical 
Method 

For the description of finite amplitude convection in rotating spherical fluid shells 
we follow the standard formulation used by many authors Zhang and Busse(1987), 
Ardes et al. (1997), Zhang(1991, 1992). It is assumed that a general static state 
exists with the temperature distribution T$ = T — f3d 2 r 2 /2 + ATr]r~ 1 (l — r])~ 2 where 
r) denotes the ratio of inner to outer radius of the shell and d is its thickness. AT 
is the temperature difference between the boundaries in the special case f3 — of 
vanishing internal heat sources. The gravity field is given by g = — ^dr where r is the 
position vector with respect to the center of the sphere and r is its length measured 
in units of d. In addition to the length d, the time d 2 jv and the temperature 
v 2 /'-/ad 4 are used as scales for the dimensionless description of the problem where 
v denotes the kinematic viscosity of the fluid and k its thermal diffusivity. The 
density is assumed to be constant except in the gravity term where its temperature 
dependence given by a = (dg/dT)/g = const, is taken into account. The basic 
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equations of motion and the heat equation for the deviation O from the static 



temperature distribution are thus given by 

d t u + u ■ Vu + rk x u = -Vn + 0r + V 2 u (la) 

V-it = (lb) 

P(d t & + u ■ V0) = (Ri + R e r]r- 3 (1 - r])- 2 )r ■ u + V 2 (lc) 



where the Rayleigh numbers Ri and R e , the Coriolis parameter r and the Prandtl 
number P are defined by 

a-f/3d G _ a-fATd 4 _ _ 2Qd 2 v 

UK UK UK 

Since the velocity field u is solenoidal the general representation in terms of poloidal 
and toroidal components can be used, 

u = V x (Vf x r) + Vw x r. 

By multiplying the (curl) 2 and the curl of equation fflal) by r we obtain two equations 
for v and w, 

[(V 2 - d t )L 2 + r^]V 2 w + tQw - L 2 = -r • V x [V x (u ■ Vu)] (2a) 

[(V 2 - d t )L 2 + rd^}w - T Qv = r-V x(u-Vu) (2b) 

where <9 t and denote the partial derivatives with respect to time t and with respect 
to the angle of a spherical system of coordinates r, 6, <fi and where the operators 
L 2 and Q are defined by 

L 2 = -r 2 V 2 + d r (r 2 d r ) (3a) 
Q = r cos ^V 2 - (L 2 + rd r ) (cos 09 r - r~ l sin (3b) 
Stress-free boundaries with fixed temperatures are most often assumed, 

v = d 2 „v = d r (w/r) = = (4) 

at r = Tj = 77/(1 — 77) and at r = r Q = (1 — rf)~ l 

but sometimes no-slip boundaries are used in which case the conditions d 2 r v = 
d r (w/r) = must be replaced by d r (rv) = w = 0. For simplicity we shall assume the 
previously used value 77 = 0.4 unless indicated otherwise. The numerical integration 
of equations (flcj) and (2) together with boundary conditions (TjJ proceeds with the 
pseudo-spectral method which is based on an expansion of all dependent variables 
in spherical harmonics for the 8, ^-dependences, i.e. 

v = ]T Vi n {r, £)P, m (cos 6) exp{im<j)} (5) 

l,m 

and analogous expressions for the other variables, w and 0. P z m denotes the associ- 
ated Legendre functions. For the r-dependence expansions in Chebychev polynomi- 
als are used. Spherical harmonics up to the order of / = 96 and up to 41 collocation 
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points in radial direction have been used. A time step as low as 10~ 6 was required 
for the computations of the small Prandtl number cases reported in section 5. For 
further details see Busse et al. (1998), Tilgner and Busse (1997). 

Benchmark comparisons for stationary solutions of equations ( TTcl) and (2) to- 
gether with conditions (jlj) when the dynamo process is included have been done in 
the past (Christensen et al, 2001). 




^ f 



Figure 1: Banana cells in a thin rotating spherical fluid shell cooled from within. The 
motion is made visible by a suspension of Bakes which become aligned with the shear. 



3 Onset of Convection in Rotating Spherical Shells 

A rough idea of the dependence of the critical Rayleigh number R ic on the parameters 
of the problem in the case of internal heating can be gained from the expressions 
derived from the annulus model (Busse, 1970) 




where r m refers to the mean radius of the fluid shell, r m = (r* + r )/2, and 9 m to 
the corresponding colatitude, 6 m = arcsin(r m (l — 77)). The azimuthal wavenumber 
of the preferred mode is denoted by m c and the corresponding angular velocity of 
the drift of the convection columns in the prograde direction is given by u c /m c . 
In figure [2] the expressions ( l6al l6cl) are compared with accurate numerical values 
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Figure 2: Critical Rayleigh number Ri c and frequency lu c (right ordinate) as a function 
of the Prandtl number P in the case rj = 0.4 for the Coriolis numbers r = 5 • 10 3 (pink), 
10 (blue), 1.5 • 10 4 (green) and 10 5 (red). The thick black line corresponds to expressions 
JSaD and (Sep . The reference value R ic (l) equals 80296, 190140, 318395, 3665919 for t = 
5 • 10 3 , 10 4 , 1.5 • 10 4 , 10 5 , respectively. For expression i6a\) -R« c (l) is given by 0.67543 • r 4 ' 3 . 

which indicate that the general trend is well represented by expressions f l6a] l6cl) . 
The same property holds for m c . For a rigorous asymptotic analysis including the 
radial dependence we refer to Jones et al. (2000). In the following we shall continue 
to restrict the attention to the case R e = and thus we shall drop the subscript % 
oiRi. 

There is a second mode of convection which becomes preferred at onset for 
sufficiently low Prandtl numbers. It is characterized by convection cells attached 
to the equatorial part of the outer boundary not unlike the "banana cells" seen 
in the narrow gap experiment of figure [TJ The equatorially attached mode actually 
represents an inertial wave modified by the effects of viscous dissipation and thermal 
buoyancy. An analytical description of this type of convection can thus be attained 
through the introduction of viscous friction and buoyancy as perturbations as has 
been done by Zhang (1994) for the case of stress-free as well as for no-slip boundaries 
(Zhang, 1995). According to Ardes et al. (1997) equatorially attached convection is 
preferred at onset for r < 17 where T\ increases in proportion to P 1 / 2 . 

In the derivation of the asymptotic relationships (6) the exact nature of the 
boundary condition for the tangential component of the velocity vector at the spher- 
ical surface does not enter. In the limit of infinite r the conditions of onset of con- 
vections are the same for no-slip and for stress-free conditions. In the following we 
shall restrict the attention, however, to the case of stress-free boundaries in order to 
avoid the complications arising from Ekman boundary layers at no-slip surfaces. 
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4 Evolution of Convection Columns at Moderate 
Prandtl Numbers 




Figure 3: Lines of constant r Jrj in the equatorial plane in the case r = 3 . 10 4 , P = 0.1 
for R = 3 • 10 5 (left) and for R = 3.5 • 10 5 (right). The four sections are a quarter of the 
vacillation period, t p = 0.0864 (left) and t p = 0.0124 (right), apart with time progressing 
in the clockwise sense. 

In general the onset of convection in rotating fluid spheres occurs supercritically. 
As long as the convection assumes the form of shape preserving travelling thermal 
Rossby waves as described by linear theory, its azimuthally averaged properties are 
time independent. In fact, as seen from a frame of reference drifting together with the 
convection columns the entire pattern is steady. A differential rotation is generated 
through the action of the Reynolds stress. The latter is caused by the spiralling cross 
section of the columns which persists as a dominant feature at moderate Prandtl 
numbers far into the turbulent regime. The plots of the streamlines r|r = const, 
in the equatorial plane shown in figure [3] give a good impression of the spiralling 
nature of the columns. 

A true time dependence of convection develops in the form of vacillations after 
a subsequent bifurcation. First the transition to amplitude vacillations occurs in 
which case just the amplitude of convection varies periodically in time as exhibited 
in the left plot of figure [31 At a somewhat higher Rayleigh number shape vacillations 
become noticeable which are characterized by periodic changes in the structure of 
the columns as shown in the right plot of figure [31 The outer part of the columns 
is stretched out, breaks off and decays. The thinning and subsequent breakup can 
actually be best seen in movie 1. The tendency towards breakup is caused by the fact 
that the local frequency of propagation varies with distance from the axis according 
to expression (6c) after 6 m has been replaced by the local colatitude 9. 

The two types of vacillations also differ significantly in their frequencies of oscil- 
lation. This is evident from the time records of the energy densities of convection 
which have been plotted in figure HI The various components of the energy densities 
are defined by 
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Movies 1, 2 and 3. 



where v refers to the azimuthally averaged component of v and v is given by v = v—v. 

As the Rayleigh number is increased further a fairly sudden transition into a 
chaotic regime occurs where convection has become strongly inhomogeneous in space 
and in time. A typical sequence of plots is shown in figure 5 which covers about 
one period of the the relaxation cycles seen in the time record for R = 3.8 • 10 5 
in figure HI In contrast to the more common relaxation oscillations encountered 
at higher Rayleigh numbers, see, for example, the time record for R = 5 • 10 5 in 
figure HI convection does not die off entirely at any time during the cycle. But 
the interaction between convection and differential rotation appears to be similar. 
As the amplitude of convection as measured by E p and E t grows the differential 
rotation generated by the Reynolds stress grows as well with just a small delay in 
time. When the differential rotation reaches a critical level the convection columns 
become disrupted and their amplitude decays. Subsequently the differential rotation 
decays as well on the time scale of viscous diffusion. It is typical for this type of 
relaxation cycles that the viscous decay is shorter than the growth time of the 
differential rotation in contrast to the relaxation oscillations at higher values of R. 
The regime of relaxation oscillations is interrupted once in a while by a return to 
the more regular regime of shape vacillations as shown in the inserted enlargement 
of the time record for R = 3.8 • 10 5 in figure HI But in contrast to the vacillation 
of the right plot figure [3] the pattern is now strongly modulated. The component 
with the azimuthal wavenumber m=l plays a dominant role in this modulation. As 
can be seen in figure 6 the pattern recurs nearly periodically in time in spite of the 
modulation. Figure 6 corresponds to the slightly higher Rayleigh number of 4 • 10 5 
where convection changes intermittently between relaxation cycles and vacillation 
oscillations. 

As R is further increased the spatio-temporal structure of convection becomes 
more irregular as can be seen in the time series for R = 4.5 • 10 5 in figure HI But 
at R = 5 • 10 5 the relaxation oscillations with only intermittent convection become 
firmly established and they continue to persist up to Rayleigh numbers of the order 
of 1 • 10 6 . They are basically the same phenomenon as found by Grote and Busse 
(2001) at P = 1 and even the period is the same, about 0.1, which corresponds 
to the viscous decay time of the differential rotation. The only difference that can 
be noted is the small precursor hump before the burst of convection. In figure 7 
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Figure 4: Time series of energy densities of convection in the case r = 3 • 10 4 , P = 0.1 
for R = 3 • 10 5 , 3.3 • 10 5 , 3.5 • 10 5 , 3.8 • 10 5 , 4 • 10 5 , 4.5 • 10 5 , 5 • 10 5 (from bottom to top). 
Red, blue and green lines indicate E™ , E{ , and respectively. The critical Rayleigh 
number for onset of convection is R c = 222518. E™ is smaller by more than an order of 
magnitude than the other energy densities and has not been plotted for this reason. 



the third plot shows convection appearing in a fairly regular pattern. Afterwards 
it decays again as shown in the fourth plot before it erupts in a chaotic fashion as 
the differential rotation has reached its minimum amplitude. Since this latter part 
of the cycle occurs on a rather short time scale, smaller intervals between the plots 
have been used for this part in figure 7. 

The main difference between convection at P = 0.1 and Prandtl numbers of 
the order unity occurs at the transition between regular and irregular patterns. A 
typical scenario is shown in figure 8. After convection has set in in the form of 
eight drifting column pairs the usual amplitude vacillations occur as the Rayleigh 
number is increased. The shape vacillations, however, exhibit a modulation with 
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Figure 5: Sequence of plots equidistant in time (from top to bottom with At = 0.05) for 
R = 3.8 • 10 5 , r = 3 • 10 4 , P = 0.1. The left column shows streamlines, = const., in the 
equatorial plane and in the middle column lines of constant radial velocity u r on the mid 
surface, r = + 0.5, are shown. The left halves of the circles of the right column show 
lines of constant ua, which is the azimuthally averaged azimuthal component of the velocity 
held. The right halves show streamlines of the axisymmetric meridional circulation. 

wavenumber m = 4 as shown in figure 9. Only every second pair of columns gets 
stretched until the outer part separates, a little earlier for the cyclonic column than 
for the anticyclonic one. Then the same process is repeated for the other pairs of 
columns such that the sequence shown in figure 9 exhibits only half a period of 
the oscillation. In addition, of course, the column pattern drifts in the prograde 
direction. Movie 2 provides a vivid impression of this type of convection. 

With increasing Rayleigh number the stretching process gets out of phase and 
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Figure 6: Modulated shape vacillations of convection for R = 4-10 5 , r = 3-10 4 , P = 0.1. 
The sequence of plots equidistant in time (At = 0.00376), starting at the upper left and 
continuing clockwise, shows streamlines, = const., in the equatorial plane. Since the 
modulation period is t p = 0.0188, the last plot closely resembles the first plot except for a 
shift in azimuth. 

an m = 1-modulation enters as shown in figure 10. The time dependence is still 
periodic, but has become more complex in that some of the separated outer parts 
become attached to the preceding column pair. To show this process in more detail 
some extra plots have been added to the time sequence of figure 10. While the 
amplitude of convection varies considerably throughout the cycle, the differential 
rotation exhibits only small oscillations as can be seen from the corresponding section 
of figure 8. 

By the time when R has reached 3.5 • 10 5 the modulated vacillations have become 
aperiodic and with increasing Rayleigh number convection becomes more and more 
chaotic. Regularity reappears only in the form relaxation oscillations as shown in 
the section for R = 7 ■ 10 5 of figure 8. Remainders of the vacillations can still be 
seen in this time record. But at R = 1 ■ 10 6 the relaxation oscillations have become 
fully established with convection occurring only in a intermittent fashion. Although 
obtained for a slightly smaller value of r, movie 3 gives a good impression of the 
onset of relaxation oscillations. 

5 Convection at Small Prandtl Numbers 

Since stellar interiors as well as metallic planetary cores are characterized by rather 
small Prandtl numbers much attention has been focused on the problem of the 
onset of convection in rotating fluid spheres in the case of small P. Zhang and 
Busse (1987) found the equatorially attached mode which is quite distinct from 
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Figure 7: Relaxation oscillations of convection for R = 5 • 10 5 , r = 3 • 10 4 , P = 0.1. The 
upper 6 plots show a time sequence (upper row left to right, lower row left to right) of 
streamlines, = const., in the equatorial plane. The corresponding lower 6 plots show 
lines of constant in the left halves of the circles and the meridional circulation in the 
right halves. The separation in time is 0.0378 for the Erst four plots and 0.0126 for the 
last three plots. 

the columnar mode discussed in the preceding section. The new mode represents an 
inertial oscillation which becomes excited when viscous dissipation is balanced by the 
energy provided by thermal buoyancy. The fact that both energies can be regarded 
as small perturbations has led Zhang (1994, 1995) to solve the problem of onset of 
convection by an asymptotic analysis. An detailed numerical study together with 
some analytical approximations can be found in the paper of Ardes et al. (1997). 
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Figure 8: Time series of energy densities of convection in the case r = 1.5- 10 4 , P = 0.5 
for R = 3 • 10 5 , 3.2 • 10 5 , 3.45 • 10 5 (from bottom to top, left side) and R = 3.5 • 10 5 , 7 • 10 5 , 
10 6 (from bottom to top, right side). Red, blue and green lines indicate E™,E{, and E{, 
respectively. The critical Rayleigh number for onset is R c = 215142. 




Figure 9: Modulated shape vacillations of convection for R = 3.2 • 10 5 , r = 1.5 • 10 4 , 
P = 0.5. The sequence of plots equidistant in time (At = 0.005J, starting at the upper 
left and continuing clockwise, shows streamlines, r& = const., in the equatorial plane. 
Since the modulation period is about 0.025, the last plot resembles the first plot except 
for a shift in azimuth. 
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Figure 10: Modulated shape vacillations of convection for R = 3.45 • 10 5 , r = 1.5 • 10 4 , 
P = 0.5. The left and middle columns show streamlines, r^r = const., in the equatorial 
plane, the right column shows lines of constant in the left halves of the circles and 
streamlines of the meridional circulation in the right halves. The left and the right columns 
show plots at the times t = (n — 1) • 0.01, for n = 1,2,3,4 (from top to bottom ), the middle 
column shows plots at the intermediate times t = (n — 1) • 0.005 

The results of these various efforts have turned out to be rather complex since the 
preferred inertial modes travel in the prograde as well in the retrograde directions 
depending on the parameters of the problem. Moreover, the azimuthal wavenumber 
m of convection does not increase monotonically with the Coriolis parameter r 
as is usually found for the columnar mode at values of P of the order unity or 
higher. In addition to the simple "single cell" inertial modes of Zhang and Busse 
(1987) multicellular modes have been found in the study of Ardes et al. (1997) 
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Figure 11: The critical Rayleigh number R c (upper plot) and the corresponding fre- 
quency uj c (lower plot) as a function of r in the case P = 10 -4 , rj = 0.3, m = 8. The 
values corresponding to the analytical expression (9) for uj c are shown with dashed lines. 



which appear to be closely related to the multicellular modes described by the Airy 
function in the analysis by Yano (1992) of the analogous problem of convection in 
the rotating cylindrical annulus. The numerical solutions for the problem obtained 
by Pino et al. (2000) also indicate the onset of multicellular convection in parts of 
the parameter space. 

Here we shall present only some examples of the variety of patterns encountered 
in a rotating sphere of a low Prandtl number fluid. In figure 11 the multiple changes 
from retrograde to prograde modes as a function of r are indicated and the corre- 
sponding changes in the patterns are shown in figure 12. To facilitate comparisons 
the wavenumber m = 8 has been kept fixed. But similar changes in the structure 
of convection are found when competing values of m are admitted. The retrograde 
and prograde modes differ little in their form. The opposite phase between toroidal 
and poloidal components of motion is the most characteristic difference as indicated 
in figure 13. The sense of outward spiralling is also opposite for prograde and ret- 
rograde modes as is evident from figure 12. But in the limit of vanishing Prandtl 
number the sense of spiralling disappears because the phase of the inertial modes 
does not vary with distance from the axis. 

Whenever r is sufficiently large the frequency omega is closely approximated by 
the frequency of the corresponding inertial modes which is given by the analytical 
expression (Zhang, 1994; Ardes et al, 1997), 

u = —^(1 ± [(1 +m(m + 2)(2m + 3)- 1 ]^). (9) 
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Figure 12: Lines of constant r0 in the equatorial plane illustrating the sequence of 
transitions for the same parameters as Fig. 11 and values of r = 5, 950, 1500, 6 • 10 5 , 
8 • 10 5 , 3.5 • 10 7 (left to right, first upper row then lower row). 
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Figure 13: Graphs in the center show the critical Rayleigh numbers R c of the competing 
prograde (green line, empty circles) and retrograde mode (blue line, Riled circles) as well 
as the actual critical value (thick red line) and the corresponding frequencies as a function 
of t in the case P = 0.001, r] = 0.2,m = 6. Contours of constant radial velocity u r (lower 
plots) and the streamlines, w = const, on the spherical surface r = 0.9 (upper plots) for 
t = 58000 and 64000 are shown on the left and right sides, respectively. 

The negative (positive) sign applies for modes drifting in the prograde (retro- 
grade) direction. Explicit values of expression (9) are given in figures 11 and 13 for 
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Figure 14: The boundary separating the retrograde (below the surface) and prograde 
mode (above the surface) for a fixed value of the wavenumber m = 6. 




Figure 15: Streamlines rj^ = const, in the equatorial plane for the case P = 0.025, 
r = 10 5 with R = 3.2 ■ 10 5 , 3.4 • 10 5 , 3.8 • 10 5 (from left to right). 

comparisons with the numerically determined values. 

For an analysis of nonlinear properties of equatorially attached convection we 
focus on the case P = 0.025 with r = 10 5 . The critical Rayleigh number for 
this case is R c = 28300 corresponding to m = 10. As R is increased beyond the 
critical value other values of m from 7 to 12 can be realized, but m = 10 and 
lower values are usually preferred. An asymptotic perfectly periodic solution with 
m — 10orm = 9 can be found only for Rayleigh numbers close to the critical value 
when computations are started from arbitrary initial conditions. On the other hand, 
perfect periodic patterns appear to be stable with respect to small disturbances over 
a more extended regime of supercritical Rayleigh numbers. Distinct transitions like 
the transition to amplitude vacillations and to structure vacillations do not seem to 
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Figure 16: Streamlines = const, in the equatorial plane for the case P = 0.025, 
r = 10 5 , R = 4 • 10 5 . Time increases by 4 • 10~ 3 from one plot to the next starting at the 
left upper corner. The eight plot is similar to the first one except for a shift in azimuth. 

exist for equatorially attached convection. Instead modulated patterns are typically 
already observed when R exceeds the critical value by 10% as can be seen in the 
plots of figure 15. These modulations are basically caused by the interaction of 
two modes with neighboring values of the azimuthal wavenumber m as indicated in 
figure 16. Here the region of small amplitude drifts in the retrograde direction with 
the difference of the drift rates for m = 7 and m = 8. The time series of energy 
densities shown in figure 17 indicate that more than two modes usually contribute to 
the dynamics of the pattern since the time dependence is not periodic as one would 
expect if only two modes interact. The computations of the time series require a 
high spatial resolution together with a small time step. The time spans indicated 
in figure 17 are sufficient for reaching a statistically steady state of the fluctuating 
components of motion since these equilibrate on the fast thermal time scale of the 
order P~ x . Only close to R c the adjustment process takes longer as can be seen in 
the case R = 3.1 ■ 10 5 where a m = 10 - pattern approaches its equilibrium state. 
The pattern corresponding to the other cases of figure 17 are the ones shown in 
figures 15, 16 and 18. The differential rotation represented by E™ relaxes on the 
viscous time scale and thus takes a long time to reach its asymptotic regime in the 
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Figure 17: Time series of energy densities of convection in the case P = 0.025, r = 10 5 , 
for i? = 3.1 • 10 5 , 3.2 • 10 5 , 3.4 • 10 5 , (from bottom to top, left side) and 3.8 • 10 5 , 4 • 10 5 , 
8 • 10 5 (from bottom to top, right side). The red, blue and green lines correspond to E™, 
e( , , respectively. E™ is measured on the left ordinate. 



examples shown in figure 17. But the differential rotation is quite weak such that 
it has a negligible effect on the other components of motion except in the case of 
the highest Rayleigh number of figure 17. Even smaller is the axisymmetric part 
of the poloidal component of motion which is not shown in the plots of figure 17. 
At higher Rayleigh numbers the equatorially attached eddies spread farther into 
the interior and in some cases become detached from the equator as can be seen 
in the plots of figure 18. In this way the convection eddies contribute to the heat 
transport from the inner boundary. But at the same time they acquire the properties 
of the convection columns which are characteristic for convection at higher Prandtl 
numbers. Accordingly the differential rotation is steeply increased at R = 10 6 and a 
tendency towards relaxation oscillation can be noticed in the upper right time series 
of figure 17. 



6 Concluding Remarks 

The results on linear and nonlinear properties of low Prandtl number convection 
in rotating spherical shells are of special interest for applications to problems of 
convection in planetary and stellar interiors. Since planetary cores often consist 
of liquid metal or of metallic hydrogen as in the case of Jupiter and Saturn, their 
Prandtl numbers are rather low. Moreover, because of high temperatures radiative 
heat transport is no longer negligible which tends to lower the Prandtl number even 
more. The latter effect dominates in stellar interiors. On the other hand, turbulence 
owing to small scale motions which remain unresolved in numerical simulations will 
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tend to equalize all effective diffusivities. But this tendency is not likely to erase 
entirely the differences in the subgrid scale transports of heat and momentum. 




Figure 18: Streamlines = const, in the equatorial plane for the case P = 0.025, 
r = 10 5 with R = 6- 10 5 , 8 • 10 5 , 10 6 (from left to right). 

Of particular interest is the effect of the Prandtl number on the dynamo action 
of convection. While the strong differential rotation associated with columnar con- 
vection tends to promote the generation of magnetic fields (see, for example, Grote 
and Busse, 2001), its near absence in the case of equatorial attached convection 
may increase the Rayleigh number required for dynamo action at a given value of 
the magnetic Prandtl number. As the Rayleigh number reaches high values more 
column like convection flows tend to appear near the inner boundary of the shell 
and are expected to promote dynamo action. Numerical studies of this problem are 
in progress and will be reported in the near future. 

The regular and modulated patterns of convection reported in this paper as 
well as the coherent structures found in the chaotic regimes of convection can in 
principle be observed in experiments through the use of techniques mentioned in 
the Introduction. Low Prandtl number fluids such as liquid metals pose severe 
experimental challenges because convection flows can not easily be visualized. On 
the other hand the temperature profiles of convection patterns can be measured 
through time records of thermistor probes implanted in the spherical boundary since 
the pattern is drifting by. See, for example, measurements of convection in mercury 
in the case of the cylindrical annulus by Azouni et al. (1986) and of convection 
in a spherical shell by Cordero and Busse (1992). For the purpose of an eventual 
comparison with experimental observation the computations reported in this paper 
should be repeated for the case of no-slip boundaries with R e ^ and Ri = 0. The 
computational effort will be more demanding but no significant qualitative changes 
are expected in that case. 
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